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Abstract 



We have observed a signal for the decay D*^ D^j at a significance 
of 4 standard deviations. From the measured branching ratio B{D*^ — > 
D+j)/B{D*+ D+TT^) = 0.055 ± 0.014 ± 0.010 we find B{D*+ D+j) = 
0.017 lb 0.004 lb 0.003, where the first uncertainty is statistical and the sec- 
ond is systematic. We also report the highest precision measurements of the 
remaining D*~^ branching fractions. 

PACS numbers: 13.20.Fc, 13.40.Hq, 14,40.Lb, 12.39.Fe 
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The decays of the excited charmed mesons, D*^ and D*", have been the subject of 
extensive theoretical [0-^ as well as experimental |]5|-|Tl| investigation. The decay of the D*^ 



via emission of a vr^ or a photon has been observed and its branching ratio well measured [|I2|. 
While the Z)*+ hadronic decays D~^n° and D*+ D°7r+) have been observed 

and are widely used to tag b quark decays, the observation of the D*^ radiative decay 
remained problematic. Both D* mesons decay electromagnetically as the result of a spin-flip 
of either the charm quark or the light quark. In the case of the D*^, the decay amplitudes 
for these two processes interfere constructively. Combined with the phase space suppression 
of the hadronic decay, this interference results in a radiative decay fraction which competes 
with the hadronic decay fraction. In the case of the D*~^, the amplitudes for the two spin-flip 
processes interfere destructively. Also, there is slightly more phase space available for the 
hadronic decay. These two conditions result in a radiative decay fraction of the D*~^ which, 
in comparison to the D*^, is significantly suppressed relative to the hadronic decay fraction. 

A great deal of interest in the radiative D*~^ decay was generated by an earlier Particle 
Data Group average of B{D*^ — > = (18 ±4)% [|14|; this value was virtually impossible 



to reconcile with theory without assuming an anomalously large magnetic moment for the 
charm quark |Q. Based on 780 pb~^ of data, a previous CLEO II analysis found an upper 



limit of 4.2% (90% CL.) for this branching fraction, a result which strongly affected not only 
the D*^ branching fractions but also many B measurements. In addition to its importance 
in measuring B meson decays, a precision determination of the D*^ branching fractions will 
provide an important test of many quark models and other theoretical approaches to heavy 
meson decays IQ . For theories built around chiral and heavy-quark symmetry (heavy hadron 
chiral perturbation theory) 0, this measurement will also provide a strong constraint on the 
two input parameters {g and /?) allowing model-independent predictions to be made on a 
wide variety of observable quantities [§]. 



The approach used in this analysis is to search in the AM^ = M{D^^) — M (D^) |15| and 
AM^ = M(Z}+7r°)-M(Z}+) distributions for D*+ events using the decay chain D*+ -D+(7 
or 7T^), K^Ti^TT^. The branching ratio 

is then determined, where N{D'^'~f)/N{D^TT'^) is the ratio of the number of D*^ de- 
cays observed in each mode, and e.,^o/e^ is the relative efficiency for finding the 7r° or 
the 7 from the corresponding D*^ decay. Assuming that the three decay modes of 
the D*+ add to unity and defining R+ = B{D*+ D°tt+)/B{D*+ D+tt°), one 
finds B{D*+ D+7) = /?+/(i?+ + Rt + 1), 13{D*+ D+tt") = 1/{R+ + i?+ + 1) and 
B{D*^ D^n^) = /{Rti + -R^ + 1). Constraints on R^ can be obtained by combining 
the known phase space for D*^ —>■ D^ti^ and D*^ D^-k^ with isospin conservation and 
the expected dependence of p-wave decay widths to yield. 




i?+ = 2 ^ =2.199 ±0.064 (2) 



(where p+o and p++ are the momenta of the and in the D*^ rest frame, respectively). 
The theoretical uncertainty in this ratio is thought to be only of order 1% H], so the error 
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is dominated by those due to the M/j* — M/j mass differences [O]. This method has the 



advantage of avoiding large systematic uncertainties due to the D meson branching fractions 
and of cancehng many systematic uncertainties associated with the reconstruction. 

The analysis was performed using data accumulated by the CLEO II detector |16] at the 
Cornell Electron Storage Ring (CESR). The CLEO II detector consists of three cylindrical 
drift chambers (immersed in a 1.5 T solenoidal magnetic field) surrounded by a time-of-flight 
system (TOP) and a Csl crystal electromagnetic (EM) calorimeter. The main drift chamber 
allows for charged particle identification via specific- ionization measurements {dE/dx) in 
addition to providing an excellent momentum measurement. The calorimeter is surrounded 
by a superconductor coil and an iron flux return, which is instrumented with muon counters. 

A total of 4.7 fb~^ of data were collected at center-of-mass energies on or near the T(4S) 
resonance. The Monte Carlo simulated events used to determine signal shapes and detection 
efficiencies were produced with a GEANT-based full detector simulation. Also, a continuum 
Monte Carlo sample (which contains roughly double the statistics of the data) was used to 
test the analysis code and methods. 

Events were required to have three or more tracks and at least 15% of the center-of- 
mass energy deposited in the calorimeter. Each of the three tracks comprising a candidate 
—>■ i^~7r+7r+ decay was required to satisfy either the K~ or tt^ hypothesis at the 2.5a 
level using dE/dx alone, and then the triplet was required to satisfy the K'n'^Ti'^ hypothesis, 
including TOP information if available, with a probability greater than 10%. The three 
tracks were then constrained to come from a common vertex, and the invariant mass of the 
triplet, under the K^ti^ti^ hypothesis, was required to be within 10 MeV/c^ (~ 1.5cr) of the 
known mass. 

Photon candidates were required to be in the best region of the calorimeter, |cos^^| < 0.71 
(where 6 is the polar angle between the EM cluster centroid and the beam axis), with a 
cluster energy of at least 30 MeV. It was further required that no charged particle track 
point within 8 cm of a crystal used in the EM cluster. If the invariant mass formed by a pair 
of photons was within 2.5a of the 7r° mass, taking into account the asymmetric tt" line shape 
and the small momentum dependence of the mass resolution, the photons were identified 
as being from a vr". The photons were then kinematically constrained to the 7r° mass to 
improve the 7r° momentum measurement. 

Photons from D*~^ —>■ decays were required to pass a lateral shower shape cut, 

which is 99% efficient for isolated photons, and not to form a tt*^ when paired with any other 
photon. The decay angle 6.y, defined as the angle of the 7 in the D*~^ rest frame with respect 
to the D*~^^s direction in the laboratory frame, was required to satisfy cos6'^ > —0.35. This 
cut helps to reduce the large combinatorial background that arises when mesons are 
combined with soft photons moving in the opposite direction. 

The combinatorial background was further reduced by requiring xd* > 0.7, where xd* 
is the fraction of the maximum possible momentum carried by the reconstructed D*~^. This 
cut also removed any contribution from B —>■ D*X events. The cuts on cos6'^ and xd* were 
determined to maximize S'^/B {S is signal and B is background) by utilizing a large sample 
of D*^ —>■ D^j events from the data as well as Monte Carlo simulated events. 

The primary difficulty in this analysis is the small size of the signal, due to the branching 
fraction, relative to a large combinatorial background and, more importantly, relative to a 
background due to D*~^ radiative decays where — > K~K~^tt~^. Unlike the D*~^, the D*'^ 
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FIG. 1. The M{K K~^tt~^) distributions for background (solid) and signal (dashed) 
Monte Carlo samples. 



- Mint) 

M{D^) mass difference is 



almost always decays radiatively. This is a major problem because the M{Dl 
mass difference is 143.97 ± 0.41 MeV ^ and the M{D*+] 
140.64 ± 0.09 MeV \T7\, so these two processes cannot be separated in the mass difference 
plot because the resolution in photon energy in the decay is ~ 6 MeV. Misidentification of 
— >■ K~K~^7i~^ as — >■ K~tt^tt^ can occur because the TOF and dE/dx information 
used for particle identification does not adequately separate K^s from tt's with momenta 
above ~ 1 GeV/c. When reconstructed under the Ktth hypothesis, the two invariant mass 
distributions partially overlap, and any attempt to estimate the fraction of under the 
peak will depend strongly on the resonant substructure of the D"^ K^K^tt~^ decay, 
as well as the momentum distribution of the -D^'s. The large contribution to the lower 
sideband further complicates the analysis by preventing the use of this sideband in a 
subtraction of combinatorial background. 

A means to veto D*'^ events, independent of the decay's resonant substructure, is to 
require that the invariant mass of the three tracks reconstructed under the K~K'^it~^ hy- 
pothesis be greater than a cut which removes all the D*~^ events. An unwanted side effect of 
vetoing D*~^ events by this method is that a cut in the KKir mass distribution greatly dis- 
torts the Ktttt mass distribution, making the relative normalization between the upper 
sideband and the signal region uncertain. Thus the use of a sideband subtraction to remove 
the combinatorial background from the mass difference plot is impossible. Fig. |I| shows the 
Monte Carlo K~K~^7t~^ mass distribution found in decays and that found in decays 
when one of the 7r"'"'s is misidentified as a K~^. Since there are two possible tracks to assign 
the mass, both combinations are tried, and the one yielding the smaller mass is plotted. 

Fig. ^(a) shows the AM^ distribution for events from the M{D~^) signal region as well as 
for those from the M{D~^) upper sideband (a region three times as wide as the signal region 
starting ^ 3a above the nominal mass). The AM^ distribution for the combinatorial 
background found in the M{D~^) sideband is quite fiat under the signal region, justifying 
the use of a first order polynomial in fitting this background. No D*~^ veto has been applied 
to the data in Fig. |3(a), so a fair fraction of the events in this "signal" are D*~^ background. 
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FIG. 2. The AM^ = M{D+j) - M{D+) distributions for (a) data before the D*+ veto has 
been applied, (b) data after the tight D*+ veto has been apphed. The large feature on the left 
of the plots is due to D*~^ — > D~^7r^ where one of the photons from the vr*' decay is not detected. 
Monte Carlo studies indicate that this decay does not contribute to the signal region. The dashed 
histograms are data taken from the upper M(D^) sideband. 



The signal was fit with a modified Gaussian, the parameters for which were obtained from a 
large Monte Carlo sample of D*^ —>■ D'^'-f events. The systematic error in the fit parameters 
was estimated by studying data versus Monte Carlo differences in the very similar decay 

Fig. ^(b) shows the AM^ signal and sideband distributions for events satisfying the D*+ 
veto requirement that M{K^ K^ti^) > 1.990 GeV/c^. Monte Carlo indicates the fraction 
of D*'^ events passing this cut is 0.002to;oo2, thus if the entire signal yield (180 ± 26 events) 
found in Fig. |^(a) were due to Z)*"*" decays, 0.4lo;4 events would be expected in Fig. 0(b). 
The fit to the AM^ distribution in Fig. |^(b) yields 68 ± 19 events. When these data are 
refit with the signal constrained to be 0.4lo;4 events, the of fit increases by 15.8, 
corresponding to a significance of 4.0 standard deviations for the D*^ D^j signal. 

The presence of D*~^ D^j decays having been established, the D*~^ veto was loosened 
to maximize S"^/ [S+B) as determined by the Monte Carlo samples. Fig. |^(a) shows the AM^ 
distribution for the events which passed the optimized D*~^ veto. The fraction of mesons 
passing the veto was determined by fitting the AM,r distribution before and after the veto 
was applied to the data. This distribution was fit with a double Gaussian plus a background 
function |18| which simulates the expected threshold behavior. Figs, ^c) and H(d) show the 
AMtt distributions, along with the fits, used to determine the D*~^ veto efficiency for 
mesons. 

The results of fitting the AM^ distribution for events which passed and for those which 
failed the D*+ veto. Figs. |(a) and |(b) respectively, were: NP""'' = 87 ± 21, A^^"*' = 95 ± 16 
(statistical errors only). Defining (A'^,) as the total number of D*^ (-D*^) in the data, 
the branching ratio was then extracted by solving the following pair of equations 
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FIG. 3. (a) AM^ distribution for data after the "optimal" M{K-K+-k+) cut (the D*+ veto) 
has been apphed. (b) AM^ distribution for the vetoed data, (c) AMjr distribution for data prior 
to the M{K-K+t:+) cut. (d) AM^ distribution for data after the M{K-K+'k+) cut is apphed. 
The dashed histograms are data taken from the upper M{D~^) sideband. 



(l-e+)iV+ + (l-e,)iV, = AT/-' 
e+A^+ + esN, = A^^-^ 

where e+ is the fraction of D"^'s which pass the veto as determined by fitting the AM^. 
distributions (A^T'' = 1650 ± 57 and A^*°*"' = 2265 ± 66, where the errors are statistical 
only), and = 0.037 ± 0.007 is the fraction of D^^^s which escape the veto as determined 
by a Monte Carlo study. We find 

-^7 ~ Nr"' -esNif'-'l ^ (4) 

= 0.055 ±0.014 ±0.010 

where the ratio of efficiencies e.,^o/e^ = 1.066 ±0.064. From this branching ratio we can then 
extract the branching fractions shown in Table |. The statistical uncertainty is dominated 
by the D^j yields, and the largest systematic uncertainty is due to variations in this yield 
when the mean and width of the signal shape was varied by an amount suggested by the 
jj*o _^ £,0^ data versus Monte Carlo comparison. A similar comparison was used to estimate 
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TABLE I. The D*(2010)^ branching fractions determined from 
the measured ratio N(D~^'y)/N(D~^TT^). The first uncertainty is sta- 
tistical, the second is experimental systematic and the third is that 
which arises from the use of Eq. (^). 



Mode 


CLEO II 


PDG 


L>+7 


(1.68 ±0.42 ±0.29 ±0.03)% 


(1-1^0:7)% 




(30.73 ± 0.13 ± 0.09 ± 0.61)% 


(30.6 ±2.5)% 




(67.59 ± 0.29 ± 0.20 ± 0.61)% 


(68.3 ± 1.4)% 



TABLE II. Estimates of the systematic 
uncertainties in the measurement of R'^. 

efficiency ratio e^o / 6% 

fitting of background 9% 

fitting of signal 13% 

veto efficiency for Df (19% on e^) 1% 

veto efficiency for D+ (2% on e+) 2% 

cos6l^ > -0.35 5% 



the uncertainty introduced by the cos6'^ cut. Table fl|hsts the various sources of systematic 
uncertainty and gives estimates for their impact on the measurement of R'^. 

In conclusion, we have observed, with 4a significance, the radiative decay of the D*~^ 
and measured B{D*+ D+-f)/B{D*+ D+n^) = 0.055 ±0.017 (statistical and systematic 
uncertainties added in quadrature). Assuming Eq. (^ and that the three branching fractions 
of the D*~^ add to unity, we find the results in Table |. The hadronic branching fractions are in 
good agreement with the current PDG averages, but with substantially reduced uncertainties 
(which are now dominated by the 3% uncertainty in -Rjr). The D*~^ radiative branching 
fraction is in good agreement with theoretical expectations and the earlier upper limits set by 



CLEO II [|ry] and ARGUS |TT|. The uncertainty in this branching fraction is due primarily 



to the large combinatorial background under the radiative signal, so one can expect that 
data taken with the new CLEO II. 5 detector, which includes a silicon tracker, to reduce this 
uncertainty further in the near future. 
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